Black-Body Radiation Correction to the Polarizability of Helium 
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The correction to the polarizability of helium due to black-body radiation is calculated near room temperature. 
A precise theoretical determination of the black-body radiation correction to the polarizability of helium is 
essential for dielectric gas thermometry and for the determination of the Boltzmann constant. We find that 
the correction, for not too high temperature, is roughly proportional to a modified hyperpolarizability (two- 
color hypeipolarizability), which is different from the ordinary hyperpolarizability of helium. Our explicit 
calculations provide a definite numerical result for the effect and indicate that the effect of black-body radiation 
can be excluded as a limiting factor for dielectric gas thermometry using helium or argon. 

PACS numbers: 51.30.+i, 06.20.F-, 47.80.Fg, 31.30.J-, 31.30.Jc, 12.20.Ds 



I. INTRODUCTION 

The recent measurement JT]] of the refractive index of he- 
lium in a microwave cavity resonator has yielded the hitherto 
most precise experimental value of the molar polarizability of 
the helium atom. 



A, 



3eo 



0.5172535(47) 



mol 



(1) 



Here, Na is Avogadro's number, eg is the vacuum permittivity 
(also called the "electric constant"), and is the static elec- 
tric dipole polarizability of helium. We also recall that is 
defined in the limit of zero density. The techniques in this ex- 
periment could lead to measurements of thermodynamic tem- 
perature or to a determination of the value of the Boltzmann 
constant. Each of these applications would take advantage of 
the fact that ad has been accurately determined theoretically 
in a series of calculations including complete leading relativis- 
tic and quantum electrodynamics (QED) a'^ corrections in 
the fine structure constant expansion ||2|S with an uncertainty 
of 0.2 ppm from the estimate of the a"* term. Here, the mag- 
nitude of the corrections is given in atomic units, i.e., rela- 
tive to the Hartree energy scale. We note that for excitation 
by low-energy radiation, as is relevant for the experiment |[l|], 
the relativistic and radiative corrections to the polarizability 
are unambiguously defined. However, for higher frequencies, 
there may be additional field-configuration dependent correc- 
tions (see Appendix E of Ref. i^). 

Black-body radiation present in the cavity resonator will 
lead to a temperature-dependent correction to the measured 
value of the helium molar polarizability, which could affect 
the interpretation of such measurements. However, the correc- 
tion due to black-body radiation was assumed to be negligible 
compared to the uncertainty of the measurement in Ref. 
(Heuristic arguments that support this assumption and indicate 
that it is also valid for argon are presented in Sec. llVlbelow.) 
The experiment |fl]] determined the polarizability of helium 
atoms through the interaction of microwaves with the atoms 
in a cavity resonator. At the same time, black-body radiation 
is present in the resonator, and it also interacts with the atoms. 



As a first approximation, the interaction of the cavity mi- 
crowaves with the atoms and the interaction of the black-body 
radiation with the atoms are independent processes and do 
not affect each other. The dominant energy shifts due to mi- 
crowave radiation on the one hand, and due to black-body ra- 
diation on the other hand, are described by the corresponding 
second-order AC Stark shifts |7 ')| and are proportional, to 
very good approximation, to the dynamic second-order dipole 
polarizability at the microwave and black-body frequencies. 
The derivation of the theoretical second-order expressions in 
both classical, time-ordered perturbation theory and in the 
field-quantized framework are contrasted against each other 
in Ref. fj^. From a quantum electrodynamic (QED) point 
of view, both photon annihilation as well as photon creation 
processes contribute to the dynamic polarizability. 

Besides the second-order polarizability, there is also a 
fourth-order effect which is due to the exchange of four 
instead of two photons with the radiation field(s). When 
the radiation is monochromatic, the total fouith-order energy 
shift is propoitional to the so-called hyperpolarizability of the 
atom lUlll . However, when the atom is simultaneously inter- 
acting with both microwave as well as black-body radiation, 
the treatment has to be modified because photon creation and 
annihilation processes of one and the same field have to be 
"matched," and this excludes some intermediate, virtual states 
of the atom+radiation field from the fourth-order expressions. 
Indeed, in generalizing the fully quantized formalism to fourth 
order, we find convenient expressions which describe the two- 
color hyperpolarizability. The resulting fourth-order energy 
shift finds a natural interpretation as a perturbation of the 
second-order dynamic energy shift due to the microwave pho- 
tons. The latter is proportional to the dynamic polarizability. 
Therefore, the fourth-order effect constitutes a correction to 
the dynamic polarizability of the atom. 

Note that some of the QED corrections to the polarizability 
of helium also involve fourth-order perturbation theory (jSnSfl, 
with black-body photon interactions being replaced by the in- 
teractions with the radiative photons. However, there is an 
important difference. E.g., for the QED coiTections to the 
Bethe logarithm, the atom only emits then absorbs virtual ra- 
diative photons, while it emits then absorbs, and absorbs then 
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emits photons with the probing electromagnetic waves. In the 
evaluation of the black-body radiation correction to the dy- 
namic polarizability, we have to take into account processes 
where the atom both emits then absorbs and absorbs then 
emits black-body a/ifi probing microwave photons. 

Our paper is organized as follows. The theoretical founda- 
tions are recalled and derived in Sec.|II] Numerical investiga- 
tions are described in Sec. [Ul] Conclusions are reserved for 
SecHYl 



II. THEORY 

In order to formulate the problem, we need to take into ac- 
count the interaction of the helium atom with two electromag- 
netic fields: (i) the microwave field used to probe the electric 
dipole polarizability, and (ii) the black-body radiation field. 
We work with a second quantized radiation field and with a 
first quantized theory for the atomic electrons. Furthermore, 
we work in the Schrodinger picture. In SI units, the Hamil- 
tonian for the helium atom coupled to an external source of 
microwaves and affected by black-body radiation is given as 

H = Ho + hojMa\j ttM + ^ hajBd^B'^B + Um + Ub- (2) 



Here, M and B are multi-indices defined as 



M = kuX 



B = knX 



B^B, 



where knj and kB are the wave vectors of the probing mi- 
crowave field and of the black-body field, and Xm.b denote 
their polarizations. We sum over the modes of the black- 
body field and assume that the microwave radiation can be de- 
scribed by a single mode. The other symbols used in Eq. (|2) 
are as follows. In the nonrecoil approximation, the helium 
atom is described by the nonrelativistic Hamiltonian 



47reo?'a 



47reori2 



(4) 



where ri2 is the distance between the electrons, e is the elec- 
tron charge, Z = 2 the nuclear charge number, eg is the vac- 
uum permittivity, and 7 (a = 1, 2) is the electron-nucleus dis- 
tance. We describe the interaction of the atom with the quan- 
tized electromagnetic fields in the length gauge. The electric 
dipole interactions of the atom with the microwave and black- 
body fields are as follows. 



2eoVA/ 

E 



?M ■ r [a\j + an 



hjjjf 



2eoVf 



CLB 



(5) 
(6) 

where the photon creation and annihilation operators are 
and a, respectively. We normalize the electric field operators 
(see Ref. ||T3|) so that the energy density of the microwave 
photon integrated over the volume Vm is equal to and 



analogously for the black-body photon. The effect of the elec- 
tromagnetic fields is assumed to be a small perturbation of Hq, 
so that a perturbative treatment of the dipole interaction of the 
atom with the electromagnetic field becomes possible. 

We first consider the second-order effect which gives the 
main energy perturbation of the helium atom due to the AC 
Stark effect. A single-mode microwave field probes the he- 
lium atom in the ground state. The ground state energy is de- 
noted as Eq and and its Schrodinger wave function is denoted 
as ipQ. First-order perturbation theory in Um gives a vanishing 
effect, and the leading correction to the energy Eq is of sec- 
ond order We then average over the polarizations and prop- 
agation directions of the microwave mode. We are interested 
in the classical limit where the number of microwave photons 
is riM 3> 1, the normalization volume is large (Vm 3> 1), but 
the ratio um /Vm remains finite, and proportional to the in- 
tensity of the microwave field. Using this formalism, one may 
easily rederive ifioll the dynamic AC Stark energy shift due to 
microwave photons. 
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^(Vo|r'7^(±c^)r^|7^o), (7) 



where 



7^(±w) 



1 



E 



Eq - EjnTf^ 



(8) 



(3) is the resolvent operator for the unperturbed helium atom, and 
= r\ +r\; that is, we are denoting the ith Cartesian compo- 
nent of the sum ri + r2 of the positions ri and of both elec- 
trons simply by r '. For repeated superscripts and subscripts, 
we use the summation convention [an example is given by the 
Cartesian superscripts i in Eq.©]. The sum over ± is nec- 
essary because we treat both photon absorption followed by 
emission as well as photon emission followed by absorption. 
The intensity Im of the microwave field in our normalization 
of the field operator is given as 



hi = 



Vm 



(9) 



For our purposes (ground state of helium perturbed by a mi- 
crowave field), we may approximate to good accuracy the 
microwave frequency by the static limit (oj = um 0) in 
Eq. (|7]i. Then, the well-known final result in second order is 
rederived. 



AE 



2eoc 



0!d : 



(10) 



where the static dipole polarizability (divided by the square 
of the elementary charge e^) in nonrelativistic limit for the 
ground state is defined by 



(11) 



Here, Ti!{Q) = l/{Ho — E^)' is the reduced Green function, 
where the reference state \^q) is excluded from the sum over 
intermediate (virtual) states. 
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We now investigate the perturbation of the dipole polariz- 
ability due to the black-body radiation. This is a fourth- 
order process in the radiation field. The atom may emit then 
absorb photons from the microwave field and also emit then 
absorb photons from the black-body field. Fourth-order per- 



turbation theory, with time-independent field operators in the 
Schrodinger picture, can then be used in order to infer the en- 
ergy shift. The result reads, when taking all combinations of 
emission and absorption into account, 



AEb = {^ 2(.V ) ^ (^ 2g V ) '^^^ TZ{ujm)Eai a\jTZ{ujB + ^m) Sb aBTl{^M)'^MaM) 
+ {~.Ba\'^{'^B)'^Ma\['Tl{uJB + ujM)'^MaM'R-{i^B)'^BaB) + {'E.Ma\jTZ{ujM)'^BaBT^{i^B + ujM)'^BaBT^{(^M)'^MaM) 

+ 2 (SBa^7^(tJB) Sbob TZ'{0) SmOm 'R-{<^m) ^uaAi) ~ {^Ba\ T^i^s) '^bo.b) (SmOm TZ^{uj)EMaM) 

- {^BO^B T^^{^b) 'B.BaB){{EMa\i T^{^m) Sj\/aM)| + ^(ajj,-WB) <^ (03,^5)^ + <^(a]^_f , -wm) ^ (oa/iWa/)^ 

+ (^{ag,-ujB) ^ {aB,i^B); (aj^^, ^ (oa/jI^m)^ , (12) 



where we use the shorthand notation 



(13) 



Furthermore, we again use multi-indices M and B as defined 
in Eq. (O. The replacement terms in Eq. (fTZt correspond to 
the exchange of black-body photon emission versus absorp- 
tion, of microwave photon emission versus absorption, and si- 
multaneous exchange of both processes. Next, we consider 
the "classical limit" of a high occupation number for both 
fields, and the low frequency limit for the microwave field, 
WA/ — i- (this approximation is always valid for microwave 
photons whose energy is low compared to the first available 
atomic transition). Furthermore, we match the summation 
over the black-body modes B with an integration over the 
frequency-dependent intensity of the black-body radiation in 
Eq. 



B 



Vb 



where Planck's law gives 



uitOBi T) dw = 



diOBu{LOB^T) 



47r^c2 cxp{hujB/kBT) - 1 ' 



(14) 



(15) 



For room temperature, the black-body spectrum has its max- 
imum at frequencies much below typical atomic transition 
frequencies. Therefore, in addition to approximating the mi- 
crowave frequency in Eq. (fTZl i by zero, we may also approxi- 
mate the black-body frequency by zero in the fourth-order po- 



larizability defined in Eq. ( fT2] i. Indeed, if we employ the ap- 
proximation ujb — > in the matrix element defined in Eq. (fT2l i 
(not in the prefactor which is proportional to w^), then we 
may even integrate over the black-body photon frequency an- 
alytically. This is explored in the following. For now, we ap- 
proach the problem by numerically integrating over the black- 
body photon frequency. The result can be written as 



AEf 



2enc 



adX = AE X: 



with the dimensionless factor 



X = 



dujB x{^B,T) . 



(16) 



(17) 



Thus, the product ad X can be viewed as an effective static 
dipole polarizability of helium in the presence of black-body 
radiation where ad is the (dipole) polarizability in the absence 
of the black-body radiation, and x is a multiplicative factor 
that gives the change in the measured value due to the radi- 
ation. The integrand xi^B,T) involves the thermal distribu- 
tion of photons. 



xi'^B,T) 



9 adeoc 



Y,g{±ujB)u{ujB,T). (18) 



The g function is obtained after summing over the photon 
modes in Eq. ( fT2] i and after taking into account all photon 
creation and annihilation processes, and reads 
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g{LOB) - - (V'ok^■7^(a.B)r'^7^(c^s)H7^'(0)r'=|V'o) - (^ok^■7^(c^B)H7^'(0)r'^7^'(0)r'^|V'o) 

- \ [{Mr'n{ujByn{ujBVn{u;By\i^G) + {Mr=n'{Qyn{uB)r''n'{oy\i^o) 

- {Mr'n{u:By\i,a) (VJok'7e'2(0)r'^>o) - {M^' 'R^ bY H^o) (^o|r-'=7e'(0)r"|7Ao)) • (19) 



The dimensionless factor x is the quantity we are looking for 
as it gives the relative perturbation x to the polarizability due 
to black-body radiation on ad according to Eq. ( fTSI l. 

To this point, we have kept SI MKSA units in all formu- 
las, as practiced by the Committee on Data for Science and 
Technology (CODATA). In calculations of atomic properties, 
it is usually more convenient to use formulas in atomic units. 
The SI MKSA static polarizability ad, the angular frequency 
lob and g{u!B) from Eq. ( fT9] l are related to their atomic unit 
counterparts ad, UJb, and ^(ws) by 



ad 

UJB 



ad 



Eh 



(20a) 
(20b) 
(20c) 



where ao is the Bohr radius and Eh is the Hartree energy. 
Then, x as defined in Eq. ( fT6l ) can be obtained as 



X = / dws x{^B,T) = -- 



± 



g{±uJB) 



ad 



cxp 



Eh _ 
Kb J 



(21) 



The atomic unit system is defined so that physical quantities 
pertaining to atoms are of order one. We can thus conclude 
that X is an effect of order a'^ which is additionally suppressed 
by the Boltzmann factor. In the following, we use atomic units 
(i.e., units with e = H = nie = 1, and where the length is 
measured in Bohr radii). 



III. NUMERICAL CALCULATION 

The crucial step in the evaluation is the calculation of the 
g function defined in Eq. (fT9] l. According to Eq. (T% . it con- 
tains both manifestly fourth-order but also products of second- 
order terms. Next, we integrate g over in the interval 
(0, oo) with a weight given by the Boltzmann distribution of 
the black-body radiation. The numerical integration is not 
completely trivial, because we have to omit poles due to reso- 
nances given by the virtual state in the denominators of the re- 
solvent TZ{ll!b)- By contrast, TZ{—ujb) has no poles. Bending 



the integration contour into the complex plane around the res- 
onances solves the problem. Moreover, all intermediate dis- 
crete states with a positive factor 1 / {Eq — E„i±ijJb) > must 
be represented very accurately as they define the position of 
the resonances. 

In practice, we are interested in temperatures that do not 
exceed the room temperature substantially (< 400 K). This 
simplifies the problem, because the weight given by the 
Boltzmann factor is exponentially suppressed on the scale of 
atomic transition frequencies. For T = 273 K, the maxi- 
mum of the black-body radiation distribution lies at wsmax ~ 
0.00244 a. u. (atomic units). If we are interested in evalu- 
ating the total effect to a relative accuracy of 1%, we may 
cut off the integration interval at ojscut = 0.05 a.u., where 
uiujBcut, 273 K) /u(wBmax, 273 K) 10"^", and the ratio is 
even smaller for lower temperatures. The function g varies 
only marginally in the frequency range relevant to the black- 
body radiation at 273K, and utBcut is still an order of mag- 
nitude less than the energy (frequency) difference of the IS* 
ground state of helium and the lowest excited states. On 
the integration path from zero to uj scut, we never approach 
the singular points in the denominators of the resolvents in 
Eq. (T% . and thus only real values of lub need to be consid- 
ered. 

The nonrelativistic wave function of the ground state V'o 
and its energy Eq in atomic units are determined based on 
the Rayleigh-Ritz variational principle. We use a basis set 
of explicitly exponentially correlated functions (see Refs. ifisl 
H and also Appendix|A| 



Ns 

E 

m— 1 



(22) 

where the parameters (a, b, c) for the ith function are 
randomly generated from an optimized box (^1,^2) x 
(_Bi , i?2 ) X (Ci , C2 ) under additional constraints Ui+bi > £ as 

well as bi+Ci > e and q 



> e, where e = ^ 2 {Ef^ - Eq) 

and Eq is the ground state energy for He+. In atomic units, e 
can be interpreted as an approximate radial momentum of the 
two-electron system that characterizes the radial exponential 
fall-off of the wave function. The minimal momentum e must 
be chosen to be large enough to be consistent with the bind- 
ing of the electrons to the helium nucleus. By requiring that 
all combinations ai + bi, bi + Ci, and ai + Ci fulfill this cri- 
terion, we ensure that the wave function falls off sufficiently 
rapidly at large ri, r2, and ri2- If a randomly generated or- 
bital fails to fulfill the requirements, we generate another one 
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until conditions are met. This method follows ideas outlined 
in Refs. ifTsi [l4ll . In order to fix ideas, we should reempha- 
size that the six boundary parameters characterizing the box, 
that is, Ai, A2, Bi, B2 Ci and C2 are subject to variational 
optimization, not the random parameters a^, bi ad 

In order to obtain a more accurate representation of the 
wave function, we use two boxes that model the short-range 
and medium-range asymptotics of the helium wave functions. 
For the calculation of the fourth-order effect which is the 
subject of this paper, matrices with a moderate number of 
2Ns = 100, 200, 300, 400 and 600 basis functions are 
fully sufficient (we use a prefactor 2 in order to clarify the — 
equal — distribution of the basis functions onto the two boxes). 

In this basis, all needed matrix elements can be represented 
as linear combinations of the integrals (see Appendix [All 



so that the dipole polarizability 



3 

2 = 1 



(26) 



can be written in terms of the dipole matrix element of the 
reference state wave function and of the perturbation \S^'p). 
Variational parameters for 6ipp are generated just as for the 
ground state, but the size Np = |iVs of the basis of P states 
is chosen to be larger than Ns used for the generation of the 
ground state. With these results at hand, it is then easy to 
calculate the other second-order element 



app = ^(V'o|r*7^2,(0)r'|V'o) = {S^p\5yJp) (27) 



r(a,6,c, ni,n2,ni2) 



'12 



-ari—br2—cri2 



(23) 



with nonnegative integers ni, 712, and ni2. Recurrence rela- 
tions for their computation are well known lITil [Tsl . The result 
for the ground-state energy extrapolated from 600 functions is 
Eo = -2.903 724 377 034118(3). The linear coefficients v, 
in Eq. ( |22] | are obtained from a solution of the generalized 
eigenvalue problem. The numerical accuracy of the results is 
estimated from the apparent numerical convergence of the ma- 
trix elements as the size Ns of the S state basis is increased. 

In view of the above considerations, we can approximate 
the black-body frequency in Eq. ( fT9] l to good approximation 
by ojb = and evaluate g{Q). This is instructive, because 
g{Q) can be broken down into distinct contributions, which 
allows us to present them separately, possibly enabling an in- 
dependent verification of the calculations if needed. We thus 
calculate first the quantity ap = {tpo\r'^TZ' {O)r'^\tpo) , which is 
directly connected to the dipole polarizability by the relation 
ad = — 2ap/3 in the nonrelativistic approximation. The re- 
solvent TZ'{0) can be replaced effectively by the sum over P 
states as in 



\4>np){4'nl 
En — En 



(24) 



where the sum, for our calculation, is only over the singlet 
states, and n is the principal quantum number For the ground 
state, all contributions from intermediate states fulfill Eq — 
En < 0. In that case, the exact representation of the P state 
component of the resolvent gives the lowest possible value for 
the polarizability, thus leading to a variational principle for the 
determination of the second-order polarizability. 

For the calculation of the g function, we also need the first 
order correction to the wave function. 



Np 



ri—bir2—Ciri2 



in^r^)], (25) 



needed for g{0). 

The two fourth-order terms that enter g{0) can be expressed 

as 



ao - {drpW\S^%) , as = {S^pyiS^s 



(28) 



where the intermediate S and D states are represented in the 
form 



\6^s)^n'sioym 



Ns 



bir2—Ciri 



in ^ r2)] , 



(29a) 



m=l 



Nd/5 



X e 



-air-i—bir2—Ciri 



pp 



(29b) 



The dominant effect due to intermediate D states comes from 
the excitation of one of the electrons, and the second case adds 
a mixture of single excitations of two electrons. For the cal- 
culation, we choose 2Ns = 2Nd = 47Vs = l^P- 

With the results for ap, app, ao and as as defined in 
Eqs. (l26T l, dZTl i and ( l28T l, (see also Table |l]i, we can proceed to 
the evaluation of 



72 



m 



(5 as + 6 ao — 3 ap app) . 



(30) 



The prefactors are determined by angular algebra Here, 
the term ap a pp is equal to the sum of the last two products 
of terms on the right-hand side of Eq. (fT9t , and apt and as 
correspond to the D and S state components of the sum of the 
first four terms on the right-hand side of Eq. ([19), in the limit 
ijJb 0. Indeed, the two-color hyperpolarizability 72 = 5(0) 
can be written in terms of just the matrix elements an, cts, 
ap and app because we set cjp ^ in Eq. ( fT9l ). 

The two-color hyperpolarizability 72 = <7(0), calculated 
here, is not to be confused with the static hyperpolarizability 7 
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TABLE I: Various quantities of interest for tiie calculation of the two-color hyperpolarizability, expressed in atomic units. The definition 
of ap, app, an and as is given in Eqs. i26\ . \2H . and l l28b . Numerical values are also indicated for the two-color hypeipolarizability 
72 = g{oJ = 0) defined in Eq. l |30t . and for the hyperpolarizability 7 defined in Eq. l |3U . 



2Ns 




(y.pp 




an 


7 


72 = g{0) 


100 
200 
300 
400 
600 


-2.074787 392 227 2 
-2.074788 259 836 5 
-2.074788 260 7312 
-2.074788 261670 8 
-2.074788 261679 1 


2.122 527 900 595 
2.122 530424 544 
2.122 530428 743 
2.122 530432 055 
2.122 530432 021 


-11.261 152615 
-11.261406459 
-11.261407 099 
-11.261407 836 
-11.261407 802 


-7.755 245 582 
-7.755 398 231 
-7.755 398 608 
-7.755 399 056 
-7.755 399 033 


43.103 075 16 

43.104 221 68 
43.104 224 56 
43.104 227 94 
43.104 227 78 


59.750 569 45 
59.752 012 04 
59.752 015 66 
59.752 019 89 
59.752 019 68 


00 


-2.074788 261682(3) 


2.122 530432 01(2) 


-11.261407 80(2) 


-7.755 399 03(2) 43.104 227 7(1) 59.752 019 7(1) 


Literature 


-2.074788 261682 (3)° 








43.104 227(1)" 





Result " was pubUshed in Ref. d], for result " see Ref. jUl- 



which is used in order to describe the fourth order perturbation 
when all four photons are from the same field fl [ItII- The 
latter can be expressed in our parameters as 

g 

7 = (5as + 6a£) - SapCkpp) . (31) 
15 

The numerical prefactors are different from the ones in 
Eq. dSOl l. Using our numerical algorithm, we may verify the 
result for the static hyperpolarizability 7 given in Ref. |Ql and 
add one more significant digit, 7 = 43.104 227 7(1) as op- 
posed to 7 = 43.104 227(1) from Ref. ij (see also TablelB. 




FIG. 1 : (Color online) xi^B, T) as a function of lcjb at T = 273 K. 
The frequency scale is in atomic units (a.u.). 

The above calculations illustrate the computational proce- 
dure for the static one-color hyperpolarizabilty 7 and its two- 
color generalization 72. For finite excitation frequency, there 
is a subtlety which deserves some extra considerations. In 
the fourth-order matrix elements in Eq. (|3CTi , the outer S state 
couples to odd parity P states by dipole transitions. The inner 
virtual states can be or D states, but they can also be even 
parity P states. Let us suppose that the reference S state has 
been coupled to a z polarized odd-parity P state. By a dipole 
transition, emitting or absorbing an y polarized photon, the 
two-electron system may couple to an even parity P virtual 
state that is proportional to the x component of the cross prod- 
uct 7^1 X r2- For — > 0, the contribution of even parity P 
states vanishes due to symmetry, but it gives a finite effect for 



100 




FIG. 2: (Color online) A plot of the dimensionless quantity x against 
the temperature T measured in Kelvin (K) is shown, with special 
emphasis on the point T = 273 K. 



nonvanishing ujb- In order to fix the notation, we define 

apE^e,jk{Srp\r'\S^P%E). (32) 

where the P even (PE) states are given by (we invoke the 
summation convention over k and I) 

I^^pe) = ^jki npE{Oy\S^p'p) (33) 

NpE 

rn—l 

We choose Np e = Ns and recall the definition of the pertur- 
bation \Sijjp) due to the ith Cartesian component of the posi- 
tion operator The definition of ap e thus involves a total of 
three resolvents of the helium atom. By numerical calculation, 
using the same incremental values for the basis sets as those 
indicated in TableHl we obtain ap^ = -0.062 951884 22(4). 
This value should not be understood as the static value of a 
dynamic forth-orderpolarizability due to even-parity P states 
(indeed, as outlined above, the contribution of even-parity P 
states to the two-color hyperpolaiizability vanishes). In our 
definition of ap e, we have arranged for the e tensors to isolate 
a nonvanishing contribution in the static limit. This approach 
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serves two purposes: (i) to implicitly define the variational pa- 
rameters used in the calculation of the resolvent TZp e in both 
the static as well as the dynamic regime, and (ii) to provide 
a reference value for the contribution of even-parity P states 
to a fourth-order matrix element that has a manifestly nonva- 
nishing static limit (but no direct physical interpretation). 

For the integration in Eq. (1211 over ujb, a 120 point Gauss- 
Legendre quadrature in the interval (OjWscut) is fully suffi- 
cient to obtain 

X = 2.742 57(1) X 10"^^ for T = 273K. (34) 

Under the approximation (/(ws) w ^(O) in Eq. ( |2TI ). which 
we can do because g varies slowly on the frequency scale 
of the Boltzmann distribution at room temperature, the in- 
tegral over ujb can by done analytically. It has a charac- 
ter similar to the Stefan-Boltzmann law with a T** temper- 
ature dependence. For clarity, we now return to SI MKSA 
units. With the atomic unit quantities ad = 1.383192 . . . and 
5(0) = 72 = 59.752019 . . . defined as in Eq. we obtain 

_4a3 2 72 [kBTV 
^ 135 ^ ad\ Eh ) 

pa 4.9372 x 10"^* ^0 . (35) 

This approximation gives x = 2.7424 x 10"^^ for T = 273 K 
with an error of order 10^"*^ compared to Eq. ( [34l i. At room 
temperature (T = 300 K), we obtain x = 3.999 x lO-^^. Fig- 
ures [T] and |2] illustrate the Boltzmann weight of the integrand 
defining x the overall dependence of x on the tempera- 
ture. As indicated in Eq. ( fTSI l, x is the multiplicative factor 
that gives the temperature-dependent effective static polariz- 
ability of helium in the presence of the black-body radiation. 

IV. CONCLUSIONS 

We have investigated the effect of black-body radiation on 
the determination of the helium molar polarizability by fourth- 
order perturbation theory in the quantized electromagnetic 
field of the probing microwave and the black-body radiation 
field. This shift of the molar polarizability is of interest be- 
cause it may affect a conceivable definition of thermodynamic 
temperature or a determination of the Boltzmann constant 
based on the measurement |[lj|. Indeed, it has been suggested 
in Ref. ifisll that the relative correction to the polarizability 
amounts to a correction as large as 2.1 x lO^^'^K^^T^ = 
16 X 10^® at T = 273 K, which would have been signifi- 
cant compared to the 9.1 x 10^^ relative uncertainty of the 
measurement iQl]. However, our fourth-order calculations in- 
dicate that the effect of black-body radiation on the measure- 
ment of the polarizability of helium reported in [ IJ] is neg- 
ligible. In fact, the result of that measurement, (A^ nioas — 

^e,thoory)Me,thcory = (1.8 ±9.1) X 10^^ SUppOrtS this COn- 

clusion. The calculation in ifisll considers only the frequency 
dependence, discussed in Ref. |fl9ll for example, of the inter- 
action of the black-body radiation with the atoms. In lowest 



order, this dependence does not affect the interaction of the 
microwave radiation with the atoms. 

Here, we have presented a fully quantized approach to the 
calculation of the black-body radiation correction to the polar- 
izability of helium, and we have evaluated all expressions nu- 
merically [see Eq. (|34l i1 as well as within a semi-analytic ap- 
proach [see Eq. (|35]|1. However, even without this formalism, 
an estimate for the effect of the black-body radiation on the 
polarizability of helium at room temperature could have been 
performed based on known literature references, as follows. 
We start from the ground-state blackbody shift 5vq = 0.11 Hz 
given by Farley and Wing ll20ll . (This value takes into ac- 
count the 300 K to 273 K temperature difference as compared 
to Ref. ll20ll .) This amounts to a relative change in the po- 
larizability of order Sv^/vYie ~ 2 x 10~^^, which affects the 
virtual states in the defining expression for the ground state 
polarizability (we denote by j/hc ~ 5 x 10^^ Hz a typical 
transition frequency in helium). Unlike the QED corrections, 
the black-body radiation effect is largest for highly-excited 
states, so a possibly larger correction would come from the 
shift of the excited-state energies due to the black-body radi- 
ation. These states enter the expression for the polarizability 
as virtual states, and the relative frequency of the virtual tran- 
sitions leads to a proportional shift of the polarizability, as 
implied by the fourth-order effect given in Eq. (fT9) . Farley 
and Wing find that the correction 5vi for each excited state i 
approaches a limiting shift of Svoo ^ 2.0 kHz as i — > oo ll20ll . 
For lower states, this gives an overestimate because the black- 
body shifts for typical excited states are larger Nevertheless, 
if all excited state energies were shifted by this amount, then 
the polarizability would undergo a relative change of order 
X{T = 273 K) = 5voo/vYic « 4 x 10"^^ which needs to 
be contrasted to our exact result (the latter is in the range of 
10~^^ and confirms the overestimation). 

In this article, we have carried out a detailed analysis of the 
effect, and obtained the precise shift of the polarizability in 
fourth-order perturbation theory, which can be expressed as 
a relative perturbation of the second order polarizability ef- 
fect proportional to a dimensionless factor x- Our analytic 
formula given in Eq. (|35T l and our numerical results obtained 
by numerically integrating the black-body spectrum confirm 
this estimate under the proviso that the shift of the lowest ex- 
cited helium states, which is of the order of a few hertz, gives 
the relevant contribution to the polarizability. In particular, 
we find a shift of x{T = 273 K) 2.7424 x lO^^* for the 
relative correction to the polarizability of the helium ground 
state. Furthermore, according to Eq. (l35t . we find a de- 
pendence for the overall shift of the polarizability with an an- 
alytic dependence of the form x ~ 4.9372 x IQ-'^^{T/KY. 
Although we have carried out numerical calculations only up 
to a temperature T ~ 400 K, an analytic estimate based on 
the known form of the black-body spectrum shows that the 
formula Eq. (|35l l should be valid to better than 1 % in a tem- 
perature range to at least 2000 K. 

Finally, we would also like to briefly discuss the size of 
the effect for other noble gases like neon and argon, which 
are of experimental interest. The two-color hyperpolarizabil- 
ity scales as {Zcsct)^^^ where Zcs is the effective nuclear 
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charge seen by the electron in the valence shell. The final re- 
sult for the shift of the polarizability due to the black-body 
radiation interaction is obtained as an integral over the dy- 
namic two-color hyperpolarizability, and its value may thus 
additionally depend on the overlap of the resonance frequen- 
cies with the maximum of the temperature-dependent black- 
body spectrum. For temperatures not exceeding 400 K, how- 
ever, the static two-color hyperpolarizability of the noble gas 
in question should provide a good estimate of the total effect, 
according to Eq. ( [35T l. We thus expect that the result for x in 
the temperature range < T < 400 K for neon and argon 
should not differ from the result given in Eq. d35t by more 
than two orders of magnitude and thus be negligible on the 
level of accuracy reached in the experiment 111 . More ac- 
curate estimates require an explicit calculation of 72 for the 
atomic reference system under investigation. 
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Appendix A: Two-electron integrals 



The two-electron integral F is defined by 



F(ni,n2,n3,a, /3,7) 



r2-7ri2 ni-1 "2-1 ns-l 



4tt / Att 



ri ''12 



(Al) 



This integral takes a very simple form when all Ui = 0, 

1 



r(0,0,0,a,/3,7) 



+ (a + 7) (/3 + 7)' 



(A2) 



The explicit form for rii > Q can be obtained by differen- 
tiation with respect to the corresponding parameter For the 
actual evaluation of F, we use compact recurrence relations 
fromRefs. ifllfisl]. 
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